Temperature Dependence of the Band Gap of Semiconducting Carbon Nanotubes 
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The temperature dependence of the band gap of semiconducting single-wall carbon nanotubes 
(SWNTs) is calculated by direct evaluation of electron-phonon couplings within a "frozen-phonon" 
scheme. An interesting diameter and chirality dependence of Eg(T) is obtained, including non- 
monotonic behavior for certain tubes and distinct "family" behavior. These results are traced to a 
strong and complex coupling between band-edge states and the lowest-energy optical phonon modes 
in SWNTs. The Eg{T) curves are modeled by an analytic function with diameter and chirality 
dependent parameters; these provide a valuable guide for systematic estimates of Eg{T) for any 
given SWNT. Magnitudes of the temperature shifts at 300 K are smaller than 12 meV and should 
not affect (n, m) assignments based on optical measurements. 

PACS numbers: 73.22.-f, 63.22.4-m, 71.38.-k 
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The temperature dependence of the band gap (Eg) is 
one of the fundamental signatures of a semiconductor, 
providing important insights into the nature and strength 
of electron-phonon (e-p) interactions. The first measure- 
ments of Eg (T) date from the dawn of the semiconductor 
era [J. Typically, Eg{T) curves shovif a monotonic de- 
crease with temperature that is non-linear at low T and 
linear at sufficiently high T 0, Q • 

Semiconducting carbon nanotubes are relatively novel 
semiconductor materials Q, with a variety of potential 
applications. Despite intensive research since their dis- 
covery, it has only recently become possible to perform 
measurements of the optical gap in individual single-wall 
carbon nanotubes (SWNTs) ll, B S i, H US- Such 
measurements, combined with information from vibra- 
tional spectroscopy, provide a route to {n, m) assignment 
of SWNTs Hi, la. Understanding Eg{T) for nanotubes 
is extremely important in this context, since experiments 
are usually performed at room-temperature and (n, m) 
assignments are often guided by comparisons between 
observed optical transition patterns and the correspond- 
ing predictions from calculations at T = K. Moreover, 
the Eg(T) signature could provide extra information for 
those assignments. 

Although it has been demonstrated that many-body 
quasiparticle and excitonic effects are crucial for the cor- 
rect description of the photoexcited states and for a quan- 
titative understanding of such optical measurements , 
the single-particle gap is still a fundamental quantity be- 
cause: (i) it is the starting point for more elaborate de- 
scriptions and (ii) trends in the single-particle gap are of- 
ten preserved by such refinements. Therefore, this work 
is devoted to describing the temperature dependence of 
the single-particle band gap of semiconducting SWNTs. 
From the calculated results for 18 different SWNTs, a 
complex dependence of Eg[T) on chirality and diameter 
is found, with an unsual non-monotonic behavior for cer- 
tain classes of tubes. This behavior arises from the differ- 



ences in sign of the e-p coupling associated to low-energy 
optical phonons. A model relation describing Eg{T) for 
any given SWNT, as a function of diameter and chirality, 
is proposed. 

The temperature dependence of Eg at constant pres- 
sure can be separated into harmonic and anharmonic con- 
tributions: {dEg/dT)p = {dEg/dT)har + {dEg/dT)anh. 

The harmonic term arises from the e-p interaction eval- 
uated at the ground-state geometry. The anharmonic 
term is due to thermal expansion. The harmonic term 
is usually the more difficult one to evaluate and it has 
challenged theorists for many years 0, 0, Q, 0, 0| . 
We follow closely the formulation of Allen, Heine and 
Cardona 1^, IB', ^ to calculate Eg (T) for semiconduct- 
ing SWNTs. In the spirit of the adiabatic approxima- 
tion we write the shift of an electronic eigenenergy 
Ai?„.k of band n and wavevector k due to static atomic 
displacements from equilibrium as a second-order Taylor 
expansion: 
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where u is a 3A^-coordinate displacement vector and D is 
the corresponding (3A^ x 3A^) Hessian matrix {N is the 
number of atoms). As usual, we express the displace- 
ments as a sum over normal modes |l7l |: 



M(a, k) = 
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ej(a, K)(a^j + a^), (2) 



where u{a, k) is the displacement along direction a of 
atom K in the unit cell, with mass M^. The electron- 
phonon couplings will be directly evaluated using the 
"frozen-phonon" scheme 18], therefore the sum only runs 
over zone-center modes j, with frequency ujj. In this ap- 
proach, a sufffciently large supercell (equivalent to a fine 
Brillouin Zone sampling) will be needed to achieve con- 
vergence. The a^j and aj are creation and destruction 



operators for phonons, and £j(a,K) are components of 
the properly normalized polarization vectors. 

We now substitute Eq. |5J| into Eq. |^ and perform 
a thermal average In the harmonic approximation, 
the linear term in u vanishes and the final result is anal- 
ogous to Ref-fTsl's: 
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where : 
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1) ^ is the Bose-Einstein occupation 



number of phonon mode j and the e-p coupling coefficient 
dEn^if^/drij is given by 
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the 18 SWNTs studied. The 
our model relation of Eq. 10, 



where Xj{a,K) = {h/AlKUJjy^^ej{a, k) are "frozen- 
phonon" displacements. From Eq. Q , this is simply the 
quadratic contribution to AEnM when the atoms are dis- 
placed along a certain frozen-phonon Xj. So, in practice, 
dEn^w/drij is calculated by performing electronic struc- 
ture calculations for ±Xj and by averaging the obtained 
energy shifts so as to eliminate the linear term. 

Structural relaxations and phonon calculations are per- 
formed using the extended Tersoff-Brenner interatomic 
potential 19]. Electronic structure is calculated using a 
multiple- neighbor, non-orthogonal tight- binding method 
|20ll2ll |. Calculations are performed for 18 different semi- 
conducting SWNT's with diameters d ~ 7.6 — 13.5 A and 
spanning the entire range of chiralities: (6,5), (7,5), (7,6), 
(8,6), (8,3), (8,4), (9,4), (9,5), (9,1), (11,1), (10,2],^(12,2), 
(10,0), (11,0), (13,0), (14,0), (16,0) and (17,0). [H 

Fig. ^shows the calculated (dots) values of AEg{T) = 
Eg{T) - Eg{0) for aU 
lines are best fits to 
to be discussed below. For clarity, results are grouped 
into four panels according to similar values of the chi- 
ral angle 9. The temperature dependence of the band 
gap is relatively small compared to bulk semiconductors: 
From the set of calculated SWNTs, the largest value of 
Eg{0) - Eg{300K) is 12 meV for the (16,0) tube, with 
idEg/dT)30QK = -5.2 x 10"^ meV/K. A strong and ap- 
parently complicated chirality and diameter dependence 
emerges: SWNTs with = {n — ni) mod 3 = 2 (full 
dots) have, in general, smaller gap shifts than i/ ~ 1 tubes 
(open dots). Most interestingly, some v ~ 2 SWNTs with 
small chiral angles show a non-monotonic gap variation 
with temperature that is positive for small T and nega- 
tive for larger T. Some of these trends, in particular the 
overall magnitude of the shifts and their j/-oscillations, 
are observed in recent photoluminescence measurements 
in suspended SWNTs 

It is challenging to explain this interesting and complex 
behavior within a unified framework. The first step is 
to analyze the contributions to Eg(T) from the different 
phonon modes. This information is contained in the e-p 
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FIG. 1: Calculated and fitted AEg{T). Full dots and lines 
correspond to = 2 SWNTs, whereas open dots and dashed 
lines correspond to v = 1 SWNTs. 



spectral function for the gap, g^F, defined as |n 



(5) 



Note that dEg/dnj, defined as the difference between 
coupling coefficients for the two band-edge states (from 
Eq. Q), can be positive or negative. 

In Fig. OJa) and (b) (lower panels) we plot g^F for 
the (10,0) and (11,0) SWNTs, which are prototype ex- 
amples of = 1 and 1^ — 2 tubes, respectively. The 
plots are restricted to the low-energy phonon branches 
(shown in the upper panels), the relevant ones to de- 
scribe Eg{T) for T < 400 K. We notice that g^F (fuU 
black line) is highly structured, reflecting the complex- 
ity of the phonon dispersion of these materials, even at 
low energies. In particular, we find that the low-energy 
optical "shape-deformation" modes (SDMs) near F pro- 
vide the most important contributions to g^F at low en- 
ergies, shown by the red dashed lines. This family of 
modes is derived, in a zone-folding scheme, from the out- 
of-plane transverse-acoustical (ZA) branch of graphene 
and they deform the circular cross-section of the tubes 
into a sequence of shapes: ellipse (the so-called "squash- 
ing mode"), triangle, square, pentagon, etc. For the 
(10,0) tube, the contribution to g^F from these modes 
is strongly negative and dominates the full g^F up to 
phonon energies equivalent to 500 K. For the (11,0), 
this contibution is positive and smaller, but still dom- 
inates g^F for T < 200 K. Beyond that temperature, 
the negative contributions from other modes start to be- 
come important. This competition explains the non- 
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FIG. 2: Dimensionless g^F (lower panels) and phonon spectra 
(upper panels) for (a) (10,0) and (b) (11,0) SWNTs. Lower 
panels: Full black lines are the total g'^F and red dashed 
lines are the contributions from the shape-deformation modes 
(SDMs) near F. The g^F is obtained from Gaussian broad- 
ening (cr — 23 K) of individual phonon contributions. Shapes 
for each SDM are also shown in the figure. Upper panels; 
Phonon frequencies (S7) measured in both cm~^ and K. SDM 
branches are highlighted in red dashed lines. 



monotonic behavior and smaller magnitudes of EgiT) 
shifts for = 2 SWNTs. 

Approximate model relations for Eg{T) are common 
in semiconductor physics Despite not being exact, 
they are useful tools for a quick assessment of Eg{T) for 
a given material. We propose a model relation for Eg (T) 
of semiconducting SWNTs in the temperature range of 
T < 400 K as a two-phonon Vina model [25| : 



AEg{T) = 



ai8i 



0202 



(6) 



where Qi and 02 (0i < 62) are effective temperatures 
for the two "average phonons" and ajQj = dEg/drij are 
their effective e-p coupling coefficients. The resulting fits 
are shown in Fig. ^as full {v = 2) and dashed {v = 1) 
lines. Note that this model is equivalent to replacing the 
highly structured spectral functions shown in Fig. |21by 
only two delta functions for each tube. Therefore, the re- 
sulting fits will necessarily be approximate. Nonetheless, 
as we see from Fig^ the largest deviation from numerical 
and model results is only 2 meV in the whole T < 400 



K range. More importantly, as we shall see below, we 
are able to obtain a completely coherent and physically 
justified dependence of the parameters aj and Qj on the 
SWNT's diameter and chirality. 

We start by considering the parameter Oi, the ef- 
fective frequency for the lowest energy phonon modes. 
These modes dictate the low-temperature (T < 100 K) 
behavior of Eg{T). These are the SDMs, derived from 
the ZA branch of graphene. In graphene, this branch 
has a quadratic dispersion at low energies and momenta. 
Therefore, within a zone-folding description, we propose 
that the energies of these modes (and therefore 61) 
should scale as the inverse square of the nanotube di- 



ameter d (in dimensionless units, d — ^/v? 



nm): 



(7) 



Therefore, Qi depends only on diameter and not on chi- 
rality. 

The remaining parameters (ai, 82 and depend 
on chirality. The task of finding this dependence is 
enormously facilitated by the following ansatz: Depen- 
dence on chirality can be expressed as polynomial ex- 
pansions /(^) of a "chirality variable" ^ = (—1)'' cos(36'). 
This ansatz is reminiscent of trigonal warping effects in 
graphite |27 ^ and is justified a posteriori by our nu- 
merical results. We consider up to 2nd-order terms: 



(8) 



where r] defines the parameter under consideration. 

The parameter ai deserves special attention because it 
must embody the intriguing sign dependence of g^F on 
1/ for the SDMs described above. We find the following 
relation: 



ai 



(9) 



The derivation of this result involves a fascinating con- 
nection between dynamical and static radial deforma- 
tions and will be presented elsewhere 28] . Similar "fam- 
ily behavior" and sign oscillations with v occur for 
static deformations induced by hydrostatic pressure p3 | 
and uniaxial stress p9l |. 

We now turn to the chirality and diameter dependence 
of a2 and 02. These parameters effectively represent 
a large number of phonon modes that start to become 
"active" at somewhat higher temperatures (between 350 
K and 500 K). In this energy range, the phonon DOS 
starts to become substantial, and we should expect that 
characteristics of particular nanotubes (chirality depen- 
dence) should gradually disappear and give rise to uni- 
versal, graphite-like features. Of course, this is only exact 
in the limit of large d, therefore, for the relatively narrow 
tubes considered here, it is wise to include some chiral- 
ity dependence (in the form of Eq. |HJ) that decays with 
increasing diameter. We therefore propose the following 
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expressions: 

92 = + :!^^^ (10) 

The overall l/d factor in a2 is simple to understand, 

since ai + a2 = limT-»oo i-S-, 0^2 is the graphite- like 
contribution to the limiting slope of Eg{T). We expect 
the high-temperature renormalization of the graphite 
bands around the Fermi point to be simply a rescaling 
of the Fermi velocity. Since the semiconducting SWNTs 
gaps are obtained, in a simplified view, by slicing the 
band structure of graphene, we expect the corresponding 
change in the gap to be proportional to the gap itself, 
i.e., to \/d. 

Equations H6I11() provide a ready-to-use recipe for es- 
timating gap shifts with temperature for any nanotube 
in this diameter range |3fl| |. The 10 parameters are ob- 
tained by best fits: A = 9.45 x 10^ K, a° = -1.70 x 10"^ 
eV/K, = 1.68 X 10-6 eV/K, 73 ' = 6.47 x lO"'^ eV/K, 
= 470 K, 7f ^ = 1.06 x 10^ K, 7!'^ = -5.94 x 10"^ 
K, B = -4.54 X lO-"* eV/K, 7^ " = -2.68 x lO'^ eV/K 
and 7^^ = -2.23 x 10"^ eV/K. 

Finally, we address the effects of thermal expansion 
in Eg{T). Available descriptions of thermal expansion 
of isolated SWNTs seem to be controversial 0, Is^ . 
We propose that, for low-dimensional structures, one 
should describe thermal expansion effects in the band gap 
in terms of anharmonic changes of internal coordinates 
(bond lengths and angles), rather than lattice constants. 
The C-C bond changes quite similarly in all carbon struc- 



tures (including diamond), although the lattice constants 
may behave quite differently ^3^]. Wc therefore use the 
well-established experimental data for the thermal ex- 
pansion of diamond '34], where the thermal expansion 
of the lattice mirrors closely the thermal expansion of 
the bond, to estimate the gap shifts due to anharmonic 
effects. This leads to very small corrections in our cal- 
culated Eg{T) values. For instance, the additional gap 
shift at 300 K for a (10,0) tube due to anharmonic effects 
is only —0.2 meV. We therefore can safely neglect these 
effects in comparison with the harmonic contributions. 

In conclusion, the calculated temperature dependence 
of the band gap of semiconducting carbon nanotubes 
shows a complex but systematic diameter and chirality 
dependence. Most gap shifts at 300 K are negative and 
small (less than 12 meV) with respect to K. Therefore, 
temperature effects should not interfere with (n, m) as- 
signments. Tubes with v — \ have generally larger shifts 
than V = 2 tubes, and some v = 2 tubes with large chi- 
ral angles even display a non-monotonic Eg{T) curve. 
All these features are explained and reproduced by a 
two-phonon model relation, with diameter- and chirality- 
depcndent parameters. 

We acknowledge useful discussions with S. Saito, A. 
Jorio, R. B. Weisman, J. Wu, P. B. Allen, J. Lefeb- 
vre, P. Schelling and Y.-K. Kwon. RBC acknowl- 
edges financial support from the John Simon Guggen- 
heim Memorial Foundation and Brazilian funding agen- 
cies CNPq, FAPERJ, Instituto de Nanociencias, FUJB- 
UFRJ and PRONEX-MCT. Work partially supported by 
NSF Grant No. DMROO-87088 and DOE Contract No. 
DE-AC03-76SF00098. 



[1] M. Becker and H. Y. Fan, Phys. Rev. 76, 1531 (1949). 

[2] O. Madelung, Semiconductors : group IV elements and 
III-V compounds (Springer- Verlag, Berlin, 1991). 

[3] M. Cardona et al., Phys. Rev. Lett. 92, 196403 (2004). 

[4] S. lijima. Nature (London) 354, 56 (1991). 

[5] A. Jorio et al., Phys. Rev. Lett. 86, 1118 (2001). 

[6] Z. M. Li et al., Phys. Rev. Lett. 87, 127401 (2001). 

[7] M. J. O'ConneU et al.. Science 297, 593 (2002). 

[8] S. M. Bachilo et al.. Science 298, 2361 (2002). 

[9] A. Hartschuh et al.. Science 301, 1354 (2003). 
[10] J. Lefebvre et al., Phys. Rev. Lett. 90, 217401 (2003). 
[11] C. D. Spataru et al., Phys. Rev. Lett. 92, 077402 (2004). 
[12] H. Y. Fan, Phys. Rev. 82, 900 (1951). 
[13] E. Antoncik, Czech. J. Phys. 5, 449 (1955). 
[14] P. B. Allen and V. Heine, J. Phys. C: Solid State Phys. 
9, 2305 (1976). 

[15] P. B. Allen and M. Cardona, Phys. Rev. B 23, 1495 
(1981). 

[16] P. B. Allen and M. Cardona, Phys. Rev. B 27, 4760 
(1983). 

[17] A. A. Maradudin, in Dynamical Properties of Solids, 
edited by G. K. Horton and A. A. Maradudin (North- 
Holland, Amsterdam, 1974), vol. 1, p. 1. 



[18] P. K. Lam and M. L. Cohen, Phys. Rev. B 25, 6139 
(1982). 

[19] D. W. Brenner et al., J. Phys.: Condens. Matter 14, 783 
(2002). 

[20] N. Hamada et al., Phys. Rev. Lett. 68, 1579 (1992). 
[21] S. Okada and S. Saito, J. Phys. Soc. Japan 64, 2100 
(1995). 

[22] We use 6-unit supercells for zig-zag tubes, which is 
enough to converge the gap shifts to within 0.1 meV in 
the K to 400 K range. Most chiral tubes already have 
a very long unit cell, therefore it is sufficient to consider 
just 1-unit cell. Also, we have verified, for the (10,0) and 
(11,0) tubes, that our empirical description of the most 
relevant (low frequency) phonon frequencies and corre- 
sponding e-p couplings compare well (within ~ 10% and 
~ 20 %, respectively) with those calculated within first- 
principles methods using the SIESTA code 35]. 

[23] J. Lefebvre et al., Phys. Rev. B 70, 045419 (2004). 

[24] R. Passler, phys. stat. sol. (b) 236, 710 (2003). 

[25] L. Vina et al., Phys. Rev. B 30, 1979 (1984). 

[26] R. Saito et al., Phys. Rev. B 61, 2981 (2000). 

[27] S. Reich and C. Thomsen, Phys. Rev. B 62, 4273 (2000). 

[28] R. B. Capaz, C. D. Spataru, P. Tangney, M. L. Cohen, 



5 



and S. G. Louie, to be published. 

[29] Y. N. Gartstcin ct al., Phys. Rev. B 68, 115415 (2003). 

[30] Wc have also calculated the temperature shifts for the 
second allowed transition {£22). The results can be un- 
derstood within the same framework as Eg (or En), pro- 
vided that the diameter dependences of ai and 02 axe 
rescaled by a factor of two and the sign factor (—1)" is 
reversed. 

[31] Y.-K. Kwon et al., Phys. Rev. Lett. 92, 015901 (2004). 



[32] P. K. SchelUng and P. Keblinski, Phys. Rev. B 68, 035425 

(2003). 

[33] N. R. Raravikar et al., Phys. Rev. B 66, 235424 (2002). 
[34] D. E. Gray, ed., American Institute of Physics Handbook, 

3rd Edition (McGraw-Hill, New York, 1972). 
[35] J. M. Soler et al., J. Phys.: Condens. Matter 14, 2745 

(2002). 



